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C/3 1 Abstract 

X. 

i-Sh , By comparing the impact of established laser smoothing techniques like Random Phase Plates (RPP) and 



Smoothing by Spectral Dispersion (SSD) to the concept of "Spike Trains of Uneven Duration and Delay" (STUD 
pulses) on the amplification of parametric instabilities in laser-produced plasmas, we show with the help of 
numerical simulations, that STUD pulses can drastically reduce instability growth by orders of magnitude. The 
simulation results, obtained with the code Harmony in a nonuniformly flowing mm-size plasma for the Stimulated 
Brillouin Scattering (SBS) instability, show that the efficiency of the STUD pulse technique is due to the fact that 
successive re-amplification in space and time of parametrically excited plasma waves inside laser hot spots is 
minimized. An overall mean fluctuation level of ion acoustic waves at low amplitude is established because of 
the frequent change of the speckle pattern in successive spikes. This level stays orders of magnitude below the 
levels of ion acoustic waves excited in hot spots of RPP and SSD laser beams. 

1 Introduction 

The concept of "Spike Trains of Uneven Duration and Delay" (STUD pulses (T)) is a novel approach that aims 
to overcome the problem strongly growing parametric instabilities in laser-produced plasmas in the context of 
inertial confinement fusion (ICF) and high energy density physics (HEDP) more generally. The method was 
conceived in 2009 by B. Afeyan [ 1| and is explained in detail in a companion article by Afeyan [2| in this 
volume. The physical implementation of STUD pulses in large laser systems offers the prospect of controlling 
the unbridled growth of electron plasma or ion acoustic waves driven by Raman or Brillouin scattering. 

So-called optical smoothing techniques, like Random Phase Plates (RPP), Smoothing by Spectral Dispersion 
(SSD), and Induced Spatial Incoherence (ISI) have been developed in the last three decades to diminish the risk 
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of hydrodynamic instabilities in fusion capsules, as well as with the hopes that residually, they may limit the 
growth of laser-plasma instabilities. The impact of these smoothing techniques on parametric instabilities, like 
stimulated Brillouin scattering, is the subject of the work presented here. Unfortunately, neither RPP nor SSD 
can efficiently inhibit growth of parametric instabilities inside potentially cooperative regimes between thousands 
of laser-intensity "hot spots" that form the typical speckle patterns of smoothed laser beams. 

"Spike Trains of Uneven Duration and Delay" require the use of phase plates such as RPP (or "CPP"), where 
the laser beam spatial profiles are a speckle pattern. However, two essential features distinguish STUD pulses 
from previous techniques (2): 

(i) instead of following a continuous profile, the laser pulse is subdivided in spikes in which the intensity 
is partially "on" and partially "off". The average intensity over the spike duration, however, equals the value 
equivalent to a continuous pulse, while the peak intensity of the spike is increased by a factor that depends on the 
duty cycle; 

(ii) the speckle patterns are scrambled in between spikes either each time or after a number of spikes. If 
the scrambling is in between each spike, they are called STUDx 1 pulses. In the case of STUDxInf (STUDxoo) 
pulses, a single RPP or speckle pattern will be deployed during the entire pulse with no speckle platter changes. 
Note that a STUDxInf pulse with 100% duty cycle is equivalent to an RPP (or CPP) beam. 

The principal idea behind STUD pulses is to control parametrically excited waves from getting re-amplified 
continuously, there where they have been excited previously. The frequent change of speckle pattern in between 
spikes has the effect that re-amplification of a plasma mode, generated in a hot spot, is very unlikely due to the 
exponential statistics of the intensities of hot spots. Our results validate the picture that the interaction volume 
tends towards being successively filled with relatively low-level ion acoustic waves (IAW). In contrast to this, 
speckle patterns in RPP beams, constantly amplify IAWs in situ, in particular in intense speckles. In case of SSD 
we can speak of repeated re-amplification which is very deleterious. 

Here we show the results of STUD pulses with a 50% duty cycle (equal "on" and "off' time intervals, on 
average, and twice higher peak intensity with respect to a continuous pulse). We also introduce a random jitter 
to the width of the on pules which is 10% in magnitude. This is in order not to excite ion acoustic waves by the 
harmonics of the periodic picket fence skeleton of the pulse. This gives rise to so-called STUD5010xl pulses. 
For details and other constructions, see |2|. We will show the advantages of STUD5010xl over the performance 
of STUD5010xInf pulses and the advantages of both of these over SSD and RPP beams. 

2 Simulations with Harmony 

We have explored the performance of STUD pulses in the case of SBS. On the acoustic spatial and temporal scales, 
it is possible to consider laser-plasma interactions in 2D that include thousands of laser speckles. This guarantees 
sufficiently good speckle statistics and can be run for upwards of a 100 ps at moderate computational expense. 
With the help of numerical simulations using the code Harmony (3), we have compared the evolution of this 
instability to cases where the driving laser is RPP, SSD, and 50% duty cycle STUD pulses with changing speckle 
patters with each spike, STUD5010xl, and with nonchanging speckle patterns, STUD5010xInf. Harmony solves 
two paraxial wave equations for the incident light field with the complex-valued amplitude E (with principal 
frequency and wavenumber wo, &o) an d for a backscattered light field amplitude E\ (with d\, ki) propagating in 
opposite direction to each other (k[ ^ -k ). In the version of Harmony used here, the ion acoustic wave (IAW) 
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relevant in the 3-wave coupling of SBS is resolved by a set of 3 (scalar) equations for continuity and for the 
momenta in z and x± of this short-wavelength component (k z ^ 2kg), (while in the version described in Ref. (3), 
a single envelope approach is used for the SBS-driven IAW) 

(<9, +u ■ V)«! +;i V ■ U! +«!V ■ u = (1), and (d, + u • V + 2v IAW )ui + c 2 s Vn, = f SBS (2), 

where no and uo stand for the average density and for advective flow, respectively, viaw for IAW damping, c s 
for the sound speed, and f S Bs for the acceleration due to the SBS coupling with both light field components, 
fsBS K "oViio^i exp(i2koz). Equations (1-2) for n\ and ui = (u[ z , ui_ x ) are solved with the help of a conservative 
Eulerian scheme, with flux correction, for the two vectors U+ = c s n\ ± n§u, \ for the (complexed-valued) density 
and velocity perturbations, n\, u\-, and u\ f . The vectors U± follow the counter-propagating characteristics of 
sound, as well as for the 3rd x-component of Eq. (2) for c s U\ iX . The IAW density perturbation n\ = (U+ + U-)/2c s 
is then used for the 3-wave coupling terms for the paraxial equations for Eq and E\ (see (3])- The usual 1st 
order equation for ti\ for 'weak coupling', when temporal growth is slower than the ion acoustic oscillations, 
corresponds to the case when the counter-propagating wave U- is not excited. Solving Eqs. (1) and (2) allow 
to access the regime of 'strong coupling' of SBS [4| that may easily appear to be relevant in high-intensity 
speckles of the smoothed beam. The dominant term for advection in the laser propagation direction in Eqs. (1-2) 
corresponds to a frequency shift Uo • V =; i2k u 0z . A flow profile is used in our simulations with ;< ;0 (z) = u(L z )[l - 
(z - L z )/u(L z )L„]. In order to seed SBS growth of possibly unstable modes with frequency ~ 2ko[c s - w,_o(z)], 
everywhere in the simulation box, a small-amplitude seed in the field E\ is introduced a the right-hand-side 
(RHS) boundary z = L z , having a broadband spectrum with 6<±> ~ 25k(,c s around the laser frequcency loq. This 
is done via a Langevin equation (d, + SojJE^L-, t) = 6<±> 5(f) with random noise generation each At, respecting 
Sa>At <k 1. The module for the 3-wave coupling, describing propagation of laser- and backscattered light, and 
for the short-wavelength IAWs driven by backscatter SBS, solved via Eqs. (l)-(2). is coupled to a nonlinear 
hydrodynamics module that describes both the plasma expansion as well as the (long-wavelength) acoustic waves 
associated with forward-SBS and self-focusing, Harmony takes into account the momentum transfer between the 
short- and long- wavelength components, see herefore Ref. [5|. The algorithm of Harmony has the options to 
take into account, or to exclude, (i) self-focusing and momentum transfer, (ii) pump depletion, and (iii) higher 
harmonics of the IAW. The simulations whose results are shown here have been carried out in two-dimensional 
(2D) geometry in a box of 4000 laser wavelengths, A , along the laser propagation direction (spatial resolution 
Az = 2Aq) and 8192/lo, in the transverse direction (Ax = Aq), with the wavelength chosen to be Aq = 0.53yum. 
For a laser / # -number of /# =8, the box length L- contains ~9 effective hot spots with an equal number of hot 
spot lengths in between them. The effective hot spot length is Lhs = nf^A® - 22QAq. The non-periodic boundary 
conditions in x result in an effective box-width of ~4000/lo (about half the entire box width), resulting in a number 
>1500 potential speckles, sufficiently numerous to avoid too strong a fluctuation from realization to realization 
(6). In all of the simulations, a laser pulse with a flat intensity profile in x, lasting 3 times the transit time of 
the box, 3L ; /c, precedes the main laser pulse. This pulse is used to initialize equivalent conditions for the noise 
in the plasma for each of the different types of main laser pulses that follow, either RPP, SSD, or STUD pulses. 
The main pulse is smoothly switched on over L-Jc. The simulated plasma at T e =2 keV, had 10% critcal density, 
Wo =0.1n c , with a inhomomeneous flow profile in the interval u 0z /c s =-3. . .-7c, with an inhomogeneity length 
(at z = L z ) of L„ = \u(L z )\L u =3000/lo- The IAW damping coefficient was (v/w)iaw =0.025. The corresponding 
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Figure 1: Reflectivity of SBS backscatter as a function of time from simulations with Harmony for different the 
different smoothing methods discussed here: RPP, SSD at 45GHz and 450GHz bandwidth, STUD5010x, and 
STUD5010 with a fixed speckle pattern STUD5010xInf. Note: upper/lower subplot in linear/logarithmic scale. 

SBS (Rosenbluth) gain for the average intensity of a RPP and SSD runs at I L =5.5x 10 14 W/cm 2 (and twice this 
value for STUD pulse peaks) in the inhomogeneous profile hence was Gsbs = 5.5 at z = L z . and the interaction 
length Li nt = 2(v/<jj)iAw[c.t + uq, z (z)]LJ uo, z (z) thus L mi = 100. . . 128Aq. For the applied STUD pulses the spike 
duration chosen was ct sp ± c =173 Aq. For the pulses following SSD a bandwidth of 45GHz was chosen. An 
alternative simulation with a 10 times larger bandwidth at 450GHz was also carried out to see the impact of faster 
smoothing, even if it lies outside any realistic modulation pace of SSD systems. 

Simulation results: The simulations with Harmony have been carried out starting from a noise level of ~ 3x 10~ 9 
of the incoming laser light pump intensity (of the late time RPP beam) in the backscattered light intensity. As 
illustrated in Fig. Q] the cases for RPP and SSD at 45GHz rapidly attain pump depletion in numerous speckles, 
attaining a SBS reflectivity around 10-20%|7|. SSD proves to show a reduction with respect to RPP on a longer 
time scale[ 8], but the efficiency of reduction is less than an order of magnitude for the realistic bandwidth value 
of 45GHz. The reflectivity for the 10 times higher bandwidth value exhibits a slower growth, but tends again to 
reflectivities above 1%. As can be seen from Fig. [2] the RPP case (and to a lower degree the SSD case) exhibits 
high-amplitude IAW structures in space due to repeated amplification in close vicinity to high-intensity speckles 
(see zones put in frames). Apparently the SSD technique, that can be associated with moving speckles within a 
limited volume, cannot avoid that IAWs are re-amplified sufficiently often by intense speckles. Compared to both 
the RPP and the SSD technique, the use of STUD pulses results in a drastic reduction of SBS backscattering. 
Figure [2] furthermore illustrates that the mean level of IAW fluctuations for the STUD pulse stays of at least 2 
orders of magnitude below the levels seen for all other cases. In particular, spatial structures that can be associated 
with the speckle pattern are washed out. It is also instructive to note that the STUD5010xInf pulse case, based 
on a fixed speckle pattern, but using the the same spikes, attains similar levels as RPP and SSD 45GHz. This 
underlines the fact that a frequent scrambling of hot spots into uncorrelated speckle patterns is an important factor 
for the success of the STUD pulse program. 

Part of the simulations have be performed on the facilities of IDRIS-CNRS, Orsay, France. Partial support 
from the programs DOE NNSA SSAA, Phase I SBIR from DOE OFES, and DOE NNSA-OFES Joint HEDP is 
acknowledged. 
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Figure 2: Snapshots in the x - z-plane, taken at t =100ps, of laser pump intensity (left column), of backscat- 
tered light intensity (center), and of IAW amplitudes (right), for the cases (top to bottom) RPP, SSD at 45GHz, 
STUD5010xInf, STUD5010xl pulses. The pump intensity is normalized to the beam average intensity (/), both 
other quantities to the level after the pulse preceding the main pulse. Note the difference in the levels (see color 
bars) for the STUD5010xl case with respect to the other cases. RPP case (upper line): two regions with particu- 
larly high IAW amplitudes (center) coinciding with intense speckle locations (left) are high-lighted in frames. 
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